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ABSTRACT 

We have conducted ultra-deep optical and deep near-infrared observations of a 
field around the z — 1.226 radio-quiet quasar 104420.8-1-055739 from the Clowes- 
Campusano LQG of 18 quasars at z 1.3 in search of associated galaxy clustering. 
Galaxies at these redshifts are distinguished by their extremely-red colours, with / — 
K > 3.75, and we find a factor ~11 overdensity of such galaxies in a 2.25 x 2.25 arcmin^ 
field centred on the quasar. In particular, we find 15-18 galaxies with colours consistent 
with being a population of passively-evolving massive ellipticals at the quasar redshift. 
They form 'fingers' in the V — K/K, I — K/K colour- magnitude plots aXV — K ~ 6.9, 
I — K ~ 4.3 comparable to the red sequences observed in other z ~ 1.2 clusters. We 
find suggestive evidence for substructure among the red sequence galaxies in the K 
image, in the form of two compact groups, 40 arcsec to the north, and 60 arcsec to the 
south-east of the quasar. An examination of the wider optical images indicates that 
this substructure is significant, and that the clustering extends to form a large-scale 
structure 2-3 /i~^Mpc across. We find evidence for a high (>50%) fraction of blue 
galaxies in this system, in the form of 15-20 'red-outlier' galaxies with I — K > 3.75 
and V — I < 2.00, which we suggest are dusty, star-forming galaxies at the quasar 
redshift. Within 30 arcsec of the quasar we find a concentration of blue (T^ — / < 1) 
galaxies in a band that bisects the two groups of red sequence galaxies. This band of 
blue galaxies is presumed to correspond to a region of enhanced star-formation. We 
explain this distribution of galaxies as the early-stages of a cluster merger which has 
triggered both the star-formation and the quasar. 

Key words: galaxies: clusters: general - galaxies:evolution - quasars:general - large- 
scale structure of the Universe 



1 INTRODUCTION 

Quasars are sparsely distributed across the universe, with 
the majority appearing to be unrelated to one another. How- 
ever quasar clustering has been observed on large scales 
in the form of Large Quasar Groups (LQGs) (e.g. Web- 



galaxies, but are see n at an earlier epoch with 0.4<z<2. It 
has been suggested (Romberg & Lukash 1994) that LQGs 



represent the progenitors of these local large-scale struc- 
tures, not only because of their comparable sizes, but also 
because of a claim that the comoving number d ensities of 



ster 1982; CrampLon, Cowley & HarLwlck 1987,1989; Clowes 



LQGs and local superclu sters are comparable (Komberg 



& Ca mpusano 1091,1994; Graham, Clowas fc Campusano 
1995), where 4-25 quasars form structures 100-200^~^Mpc 
across. They are thus comparable in size to the largest struc- 
tures seen at the present epoch, such as the 'Great Wall' of 



Kravtsov & Lukash 1996) 



Observations of the galaxy environments of individual 
quasars lends credence to the hypothesis that LQGs trace 
superclusters. Quasars at z ~ 0.5 are known to favour young, 
low to moderately rich clusters, with richness increasing with 
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both redshift and radio-loudness (EUingson, Yee & Green 
1991). At hig her redshifts, a sur vey of 31 1<2:<2 radio- 



loud quasars (Hall & Green 1998) indicates that they are 
on average found in moderately rich {RAbeii ~ 1-5) clus- 
ters. Optical and narrow-band (O ll) observations of fields 
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centre d on 11 z ~ 1.1 quasars from the Crampton et al. 
(1989) LQG show an excess of galaxies around all but one 
of the quasars (Hutchings, Crampton & Johnson 1995, here- 
after HCJ; Hutchings, Crampton & Persram 1993, hereafter 
HCP). In particular they find excesses of blue and emission- 
line galaxies around several of the quasars indicating re- 
gions of significant star- formation, a phenomenon probably 
related to the fact that both quasars and star-formation re- 
quire the disruption of the galaxy causing a collapse of sig- 
nificant amounts of gas. They find no evidence for a coher- 
ent structure connecting the quasars, but star-formation is 
likely to be localised, rather than occurring simultaneously 
across large scales, as a cluster mass halo collapses or during 
a cluster merger. 

To determine whether LQGs trace superclusters re- 
quires observations capable of identifying the quiescent 
galaxies that marli out any associated supercluster, in par- 
ticular the massive ellipticals which dominate cluster cores. 
These are both the most luminous and the reddest galaxies 
in nearby clusters, and form a homogeneous population with 
very tight line ar colour-magnitude (C-M) r elations known as 
red sequences (Bower, Lucey & Ellis 1992), indicative of old 



(12-13 Gyr) stellar populations. This red sequence has been 
followed in the optical for clusters out to z ~ 0.9 (Aragon- 
Salamanca et al. 1993; Stanford, Eisenhardt & Dickinson 
1998) and is consistent with the passive evolution of galax- 
ies that formed in a monolithic collapse at 2>3 ( Eggen et 
al. 196^. At 2>1 these galaxies should be characterised 
by extremely-red optical-NIR colours (J — iC' ~ 4) as the 
4000A break is redshifted into the I band. Good contrasts 
of z>l clusters over the background are possible as field 
7 — ~ 4 galaxies appear rare: the Hawaii K-band survey 
finds no galaxies with J — A" > 4 for < 18 over an area of 
86.7arcmin^ (Cowie et al. 1994; Songaila et al. 1994). 

Several z ~ 1.2 clusters have been found by searching 
for galaxies with I — K ~ 4 in fields around both targeted 
high-redshift AGN (Dickinson 1995, hereafter D95; Yamada 
et al. 1997; Tanaka et al. 2000a, hereafter TOOa) and regions 
of extended X-ray emission (Stanford et al. 1997, hereafter 
S97; Rosati et al. 1999, hereafter R99). These galaxies are 
found to have the optical-NIR colours {I—K ~ 4, R—K ~ 6) 
expected for passively-evolving galaxies which are 2-3 Gyr 
old, in good agreement with the monolithic collapse model 
predictions. D95 also finds a red sequence for the z = 1.2 
cluster around the radio-galaxy 3C324 with R — K ~ 5.9 
and an rms scatter of only 0.07 mag, suggesting that these 
galaxies formed within 300 Myr of one another. 

In this study we aim to establish the galaxy environ- 
ment of LQGs in a manner unbiased with respect to galaxy 
type. We have targeted a 30 x 30 arcmin^ field containing 
three quasars from the z ~ 1.3 Clowes & Campusano LQG 
(Fig. ^) for ultra-deep optical observations complemented 
with deep NIR observations of selected subfields, capable of 
detecting the passively-evolving galaxies that should mark 
out any coherent large-scale structures. We present here the 
first results from this study, using ultra-deep V and / im- 
ages and a K image centred on the z=1.226 radio-quiet (not 
detected at the ImJy level by the FIRST VLA 20cm sur- 
vey) quasar 104420.8-1-055739. We have detected clustering 
of galaxies with the extremely-red colours expected of quies- 
cent ellipticals at the quasar redshift. Throughout the paper 
we adopt go ~ 0.5 and Ho = lOOftkms"^ Mpc~^. 
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Right Ascension (J2a00) 

Figure 1. The Clowes-Campusano LQG with each quasar la- 
belled by its redshift. The 30 X 30 arcmin^ region targeted for 
BTC imaging is indicated by the box, and the quasar for which we 
have K imaging is circled. The boundaries of the plot match the 
boundaries of the AQD survey of ESO/SERC field 927 (Clowes 
& Campusano 1991,1994; Clowes, Campusano & Graham 1999). 
Note that this plot and Figures |5| ^ and ^ have east to the right. 



2 OBSERVATIONS 

The 30 X 30 arcmin^ field was observed through V and I fil- 
ters using the Big Throughput Camera (BTC) on the 4-m 
Blanco telescope at the Cerro Tololo Inter- American Obser- 
vatory on April 21/22 and 22/23 1998. The BTC is made up 
of 4 2048 X 2048 CCDs which have pixels of size 0.43 arcsec 
giving a field of view for each CCD of 14.7 x 14.7arcmin^. 
The CCDs are arranged in a 2 by 2 grid and are separated by 
gaps of 5.6 arcmin. To obtain a contiguous image it was nec- 
essary to shift the telescope between exposures. 16 exposures 
arranged in a 4 by 4 grid, with each pointing offset by the 
gap width to adjacent pointings, were found to give a con- 
tiguous and almost uniform coverage over a 30 x 30 arcmin'^ 
field, with 90% of the field covered by exactly 9 of the 
16 exposures. The images had large relative offsets, and so 
that they could be registered, the effect of the distortions 
(produced by the camera optics) were removed using models 
based on observations of astrometric fields. These distortions 
are significant due to the wide field of the camera, and were 
of the order of 60 pixels at the corners of the CCD array. 
After the removal of distortions, the rms errors of registra- 
tion between images were typically ~ 0.05 arcsec. Details of 
the observing and reduction procedures for the whole BTC 
image will be described elsewhere. Most of the reduction 
processes were performed in the usual way using iraf tools. 

K imaging was obtained for a 2.25 x 2.25 arcmin'^ 
field centred on the radio-quiet quasar at 10*^ 46™ 56.70", 
-f05°41'50.5"(J2000) using the UFTl camera on the 3.8-m 
UKIRT telescope in March 1999. This has a field of 
1.5 X 1.5 arcmin^ with a pixel size of 0.09 arcsec. The field 
around this quasar was observed in preference to the other 
two because of its more populous immediate environment 
(within 5 arcsec), suggesting that it was the most likely to 
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Filter 
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Seeing 


Completeness 
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time 


FWHM 
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(subficld) 


(sec) 


(arcsec) 


90% 50% 


K (centre) 


3780 


0.60 


19.25 19.95 


K (edge) 


1260 


0.60 


18.82 19.53 


/ 


16200 


1.15 


24.68 25.39 


V 


8100 


1.30 


25.40 25.95 



Observed 



Expected 



Table 1. Photometric properties of the images, including expo- 
sure times and completeness levels. The V, / figures are for the 
area covered by the K image, rather than for the whole BTC 
field. 



be in an interacting system. No cluster-scale environmental 
factors were considered, and so, for the overall clustering 
statistics, the selection can be reasonably regarded as ran- 
dom. A standard 9 point jitter pattern was used, with each 1 
minute exposure offset by 20 arcsec relative to adjacent ex- 
posures. A dark frame was obtained between each set of jit- 
ters, and subtracted before self-flattening the set. Although 
the 20 arcsec wide strips at the edges of the image have only 
one third of the exposure time, they have been retained with 
each source checked visually, and separate magnitude limits 
determined. 

The combined K image was convolved to the same see- 
ing as the / image, so that galaxy colours could be deter- 
mined using a single fixed aperture, and then registered with 
the optical images. Photometric calibration of the V and / 
images onto the Landolt system was obtained using Lan- 
dolt standard stars at varying airmasses, and a UKIRT faint 
standard was used for the K image. 

Object detection wa s carried out using SExtractor 



(Berlin fc Arnouts 1996) for objects with ~ FWHM^ con- 



tiguous pixels over the la detection threshold. Completeness 
levels for each filter (for the region covered by the K image 
rather than the whole BTC field) are given in Table ^. They 
were determined by dimming a bright galaxy to a specified 
magnitude and adding 100 copies to the image at random 
positions, and then processing through SExtractor in the 
usual way. The 50% and 90% completeness levels are esti- 
mated as the magnitude for which 50 and 90 galaxies out of 
the 100 are recovered. The total magnitudes were taken to 
be the MAGJ3EST output from SExtractor, and colours 
were determined using fixed apertures of diameter 2.5 arc- 
sec. 

Star-galaxy separation of each source was performed us- 
ing a combination of morphology (using the stellarity classi- 
fier in SExtractor) and a comparison oi V — I and I — K 
colours with model star and galaxy tracks. The stellarity 
classifier uses a neural net trained with a set of artificial 
stars and galaxies to produce a 'probability' that a source is 
stellar. Stellarities of sources in both V and I images were 
obtained and the weighted mean used for classification. Most 
of the sources classified as stars morphologically, also had 
colours that lay near the model star tracks, but a number of 
sources had the colours of blue galaxies, and were reclassified 
as such. 



Magnitude 


Total 


I-K>3.75 


Total 


I-K>3.75 


15 < /sT < 16 


1 


1 


0.533 





16 < K < 17 


3 


1 


2.186 





17 < K < 18 


13 


5 


4.204 


0.050 


18 < < 19 


23 


16 


11.36 


1.939 


19 < < 20 


39 


23 


21.81 


1.896 


Total 


79 


46 


41 


4 



Table 2. A comparison of total and extremely-red (I — K > 
3.75) galaxy counts in the K image, binned by magnitude, and 
those expected for the same-sized region in the field (Songaila 
et al. 1994). The numbers include those observed only in the K 
image, with the 2 companions to the quasar assumed to have 
I ~ K < 3.75, and the other three sources to have I — K > 3.75. 



3 RESULTS 

3.1 Galaxy Counts 

In total 100 visually-verified sources were detected in the K 
image, of which 95 had counterparts in the / image within 1 
arcsec, and 79 had counterparts in both V and I images. 14 
of the 79 sources were classed as stars, leaving 86 galaxies 
in the 2.25 x 2.25arcmin^ field. Of the 5 sources detected in 
only the K image, 2 are the close companion galaxies of the 
quasar which could not be separated from the quasar point- 
spread function by SExtractor in the V and / images, 
and the 3 others must have / — > 3.75 to be undetected 
in /. 

Table ^ shows the comparison of galaxy counts seen 
i n the K image with those expected from a field region 
( ^ongaila et al. 1994 ). Excesses of galaxies are seen in all 
magnitude bins, even in the last bin where incompleteness 
should reduce any excess. If we limit ourselves to if < 19 
(an L*if galaxy at z — 1.226 has K ~ 19.05) where we can 
be reasonably confident of both completeness and photom- 
etry, then a 3.5a excess is observed (accounting for the ef- 
fect on galaxy statistics of the two -point angular cor rela- 
tion function, uj{e) = 1.136l(")~° '* ( [Roche et al. 1999| ), for 
if < 19 galaxies), with 40 galaxies observed whereas only 
19 would be expected. By considering only those galaxies 
with / — if > 3.75 then the excess is much more significant, 
with 23 galaxies observed instead of the expected 2. Even 
though the clustering amplitude of extremely-red galaxies 
i s a factor ten grea ter than that of if-selected field galaxies 
(Daddi et al. 2000), we still observe a 9a excess. It is clear 
that the total excess is due entirely to these extremely-red 
galaxies, indicating a likely cluster at z>0.8, and we consider 
all / — if > 3.75 galaxies to be cluster members. 



3.2 Galaxy Colours and Photometric Redshift 
Estimates 

In Fig. ^ we show the I — K versus V — I diagram for galax- 
ies in the if image. To estimate the photometric redshifts 
of each galaxy from its VIK colours, and to produce the 
model colour tracks of Fig. |[ we have used the hyperz 
code of Bolzonella, Miralles & Pello ( 200C ) , which uses the 
Bruzual & Chariot evolutionary code (GISSEL98, Bruzual 
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Figure 2. I — K against V — I colour-colour diagram of all galaxies in the K image. Those galaxies whose colours are well described by 
the burst model at 2 ~ 1.2 are shown as filled squares. The four galaxies whose redshift probability distributions are shown in Fig. ^ are 
indicated by superimposed, labelled diamonds. For comparison, model tracks are shown for an instantaneous burst at zjarmation = 4.5 
(solid line), and exponentially-decaying star-formation rate models with time scales of 1 Gyr (dashed line) and 5 Gyr (dot-dashed line). 
Each track shows the colour evolution from z=1.5 to z=0 with crosses at 0.1 redshift intervals. The effect of internal extinction at 
z=1.226 is shown as an arrow in the bottom-right corner. The dotted lines indicate the selection criteria used within this paper: the 
cluster candidates are selected to have I — K > 3.75; with further separations into red {V — I > 2.00) and red outlier {V — I < 2.00) 
galaxies; blue galaxies are identified by V — / < 1.00. 



& Chariot 1993) to build synthetic template galaxies. It has 
stellar populations with 8 star-formation histories, roughly 
matching the observed properties of local galaxies from E 
to Im type: an instantaneous burst; six exponentially de- 
caying SFRs with time-scales r from 1 to 30 Gyr; and a 
constant star-forming rate. The m odels assume solar metal- 
licity and a Miller & Scalo IMF ( 1979 ), with intern al red- 
dening considered through the Calzetti et al. (MOO) model 
with Ay allowed to vary between and 1 mag. The hyperz 



software then produces a photometric redshift probability 
distribution through a chi-squared minimization process, al- 
lowing for all possible galaxy ages, star-formation histories 
and Avs. This approach of determining a range of compat- 
ible redshifts, rather than a single best-fitting redshift for 
a galaxy, is more appropriate in this case given the limited 
colour information available. 

The model curves of Fig. ^ correspond to stellar pop- 
ulations formed in an instantaneous burst ai z = 4.5 (solid 
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Figure 3. Colour-magnitude diagrams of galaxies in the K image. 
The solid symbols represent those galaxies detected in all three 
bands. The empty symbols represent those galaxies detected in 
I and K only. Galaxies whose colours are well described by the 
burst model at ^ ~ 1.2 have square symbols. The four galaxies 
whose redshift probability distributions are shown in Fig. Q arc 
indicated by diamonds. The dash-dotted lines indicate the 50% 
completeness levels for each filter in the centre of the image, with 
the dotted line indicating the completeness level in K for the edge 
of the image. 



line), and stellar populations with exponentially-decaying 
SFRs with time-scales (r) of 1 Gyr (dashed line) and 5 Gyr 
(dot-dashed line), and are thought to approximate the colour 
evolution of massive elliptical and disk-dominated galaxies. 
The galaxy ages for each model are calculated for an ^ = 0.5, 
qo = 0.5 universe. The I — K colours of the model tracks 
increase monotonically with redshift to z ~ 1.3, and it is 
clear that 7 — A" ~ 4 is a good indicator of galaxies at 
with predominantly old stellar populations. The divergence 
of the disk and burst model tracks at high redshifts is a clear 
indication of the effect of recent star- formation on the 1/ — J 
colour in particular. 

In Figs ||(a)-(d) we show the redshift probability dis- 
tributiqns, P(z), for 4 galaxies which are compatible with 
being 



L the quasai ledshift, incompatible with being at 
2:<0.8, and which represent the differing classes of galaxies 
discussed in the text below. 



Red Sequence Galaxies 

There are 15-18 galaxies (hereafter labelled 'red') whose 
extremely-red colours {I — K > 3.75, V — I > 2.00) are 
consistent with being passively-evolving galaxies at z ~ 1.2 
(indicated by squares in Fig. These red galaxies can be 
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Figure 4. Redshift probability distributions for 4 galaxies in the 
K image based on their VIK colours. The best fitting model is 
described by its redshift zjes^ and star-formation history, which is 
either an instantaneous burst or an exponentially decaying star- 
formation rate with time-scale t. The vertical dashed line indi- 
cates the redshift of the quasar. Distributions (a) and (b) cor- 
respond to two galaxies which help make up the observed red 
sequences in both V — K/K and / — K/K plots. Galaxy (c) is a 
'red outlier' with a similarly red I — K colour to (a) and (b), but 
a much bluer V — I colour, indicating more recent star-formation. 
Galaxy (d) is one of the galaxies from the blue 'band' and, due 
to its red I — K colour, is best fit by a 2^1 disk galaxy. Note that 
the central dips in the probability distributions (c) and (d) are 
only artifacts due to the limited number of star-formation regimes 
available. 



seen (as squares again) in the V — K and I — K against K 
C-M diagrams of Fig. ^ to form 'fingers' at 1/ — ii" ~ 6.9 
and I — K ~ 4.3, comparable in form to the red sequences 
observed in lower redshift clusters, and in colour to other 
clusters at the same redshift. The z ~ 1.2 clusters of R99, 
S97, and D95 are all observed to have red sequences at 
R~ K ^ 5.9, each with ^4 members spectroscopically con- 
firmed as being at the cluster redshift. The red sequence 
galaxy colours can be compared directly for the R99 cluster 
where I — K colour data exists, and the 4 spectroscopically 
confirmed red sequence members all have 4 < I — K < 4.4, in 
good agreement with ours. The colours of these red sequence 
galaxies are well-fitted by the passively-evolving monolithic- 
collapse models of elliptical galaxies, and both their s pectra 



and morph ologies are similar to present-day ellipticals (Dick- 



inson 1997). Red sequence galaxies are usually the most lu- 
minous cluster members, the massive ellipticals, and this 
also appears to be the case here, with the 10 brightest galax- 
ies with I — K > 3.75 also having V — K ^ 6.9. Given that 
these are both the reddest and most luminous galaxies, these 
should provide the tightest redshift estimates, and Fig. ^ 
shows that the third brightest of these galaxies (G3) has 
1.2<z<1.4. 

The mean I—K colour of the red galaxies is 4.25 with an 
intrinsic dispersion of 0.15 mag, and is best fit by a 2 Gyr old 
burst with a corresponding age dispersion of 400 Myr. The 
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Figure 5. Spatial distribution of galaxies in the field of quasar 104420.8+055739. The different symbols indicate the galaxies' colour, 
with circles (squares) indicating those galaxies detected (not detected) in I. Solid symbols indicate the red {I — K > 3.75, V — I > 2.00) 
galaxies which could be early-type galaxies at the quasar redshift. Half-filled symbols indicate the red outlier {I — K > 3.75, V — I < 2.00) 
galaxies which are likely to be high redshift galaxies with some recent star-formation. The star symbols indicate the blue (V ~ I < 1.00) 
galaxies which are probably undergoing significant star-formation. The size of the symbols (except stars) indicate the Kt^tal magnitude. 



colour distribution seen in the red sequence is comparable 
to that seen in the cluster of TOO (0.22 mag in J? — K ), but 
is larger than that seen in the z ~ 1.2 3C 324 cluster (D95) 
which has an rms scatter of 0.07 mag in R — K, suggesting 
that our cluster is less dynamically evolved than that of 3C 
324. 



Red Outlier Galaxies 

Only half of the excess of I ~ K > 3.75 galaxies is accounted 
for by the red sequence members, and we find comparable 
numbers (15-20) of galaxies (hereafter labelled 'red 

outliers') with both I~K > 3.75 and V-I < 2.00, which ap- 



pear to fit neither the passively-evolving nor exponentially- 
decaying star-formation galaxy models (Fig. |^ . 

Galaxies with similar colours and magnitudes have been 
observed in other z ~ 1.2 clusters (e.g. TOO; Kajisawa et al. 
1999; 2000). They have also been observed at 20<ii:<22 
in deep optical/NIR surveys (e.g. Moustakas et al. 1997), 
and appear common (several per square arcminute) at these 
fainter magnitudes. Much of the discussion (see e.g. Mous- 
takas et al. 1997) of these objects has been limited to specu- 
lation due to lack of spectroscopic observations, but a widely 
held view is that they are probably high-redshift objects 
{l<,z<^2) which are undergoing significant star- formation 
and whose extremely-red I — K colours are caused by a 
combination of dust and dominant old stellar populations. 
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Figure 6. Estimated density distribution of tiie red galaxies in 
the field of quasar 104420.8+055739. The first contour corre- 
sponds to a density of 1 galaxy per square arcminute, and the 
separation of successive contours is also 1 galaxy per square ar- 
cminute. The quasar and red galaxies are marked as in Fig. ^ 
and the labels correspond to the galaxy IDs of Table ^, which arc 
numbered in order of increasing K magnitude. 



There is also some spectroscopic evidence that some z ~ 1.2 
cluster ellipticals are undergoing star-formation, such as O ll 
emission-lines seen in object ^4 of R99 and object #237 
of S97, both of which appear bluer in optical colours than 
the other red sequence galaxies in the clusters. Some ex- 
treme members of this population have been observed with 
I ~K > 6 (e.g. Hu & Ridgway 1994), and one has since been 
spectroscopically confirmed as an ultraluminous infrared 
galaxy at z=1.44 with significant ongoing star-formation 
that is heavily obscured by dust (Graham & Dey 1996; Dey 
et al. 1999). Given that both star- formation and dust are 
likely to have affected the colours of these 'red outliers' sig- 
nificantly, it is not possible to constrain the galaxy redshifts 
beyond z^l (Fig. kr), but their prevalence in the vicinity of 
other high-redshift clusters, and the relative rarity of field 
I — K > 3.75 galaxies, suggests that many are associated 
with the cluster. 



3.3 Spatial Distribution of Galaxies 

The spatial distribution of galaxies detected in the K image 
is shown in Fig. |^, with solid symbols used to differentiate 
the red galaxies from the remainder. 



Red Sequence Galaxies 

The red galaxies appear to be distributed across the K image 
(1.29/i~^Mpc at z=1.226) with no concentration towards the 
quasar. However, the ten most luminous red galaxies (which 
we indicated earlier are probably massive ellipticals at the 
quasar redshift) are concentrated in two compact groups, 
one towards the top-centre of the image (along with a num- 
ber of fainter members), and the other in the south-eastern 



corner. This suggests that the galaxy excesses and red se- 
quences could be due to two clusters at similar redshifts. 
Assuming that the brightest of the red galaxies, Gl, is a qui- 
escent galaxy at the quasar redshift, then it has L « 14L|(-, 
and is more luminous in K, by almost a magnitude, than 
any of the galaxies from other z ~ 1.2 clusters (D95; S97; 
R99). It is also a radio emitter, being the only source in the 
K image detected by the VLA FIRST 20cm survey, having 
an integrated flux of 2.63 ± 0.15mJy, which suggests that it 
has an active nucleus. It is common for the brightest cluster 
galaxy (which is what we assume this is) to also be a radio 
source, although it is not as spectacular an example as 3C 
324 which is at a similar redshift, but is 1000 times more lu- 
minous in the radio. If this galaxy is comparable to 3C 324 
then it may also display narrow-line emission, in particular 
On and Mgll, and this may explain why it appears bluer 
than the other red sequence galaxies, as the emission-lines 
boost the optical flux. The brightest cluster galaxy is usu- 
ally located near the cluster centre, and it appears to be the 
case here too if Figs. ^ and ^ are compared. 

To examine the significance of any substructure for the 
red galaxies the non-parametric cluster analysis method de- 
scribed by Pisani (1993,1996) is applied. This is based upon 
estimating the underlying probability density, F{x), of a 
data sample Dm = {xi,X2, ■■■,xn} using an iterative and 
adaptive kernel method. Each data point, Xi, is initially 
smoothed by a Gaussian kernel of width a to produce a 
pilot fixed-kernel estimate of the local density, fp{xi), at 
each point. The optimal value for the smoothing parameter 
a is determined by minimizing the integrated square error, 
ISE{fp), of the estimation of F{x). It is possible to show 
(see Silverman 1986) that minimizing the cross-validation, 
M{fp), a function related to ISE{fp), and that can be writ- 
ten as a sum of functions dependent only on the Xi, is equiva- 
lent to minimizing ISE{fp). The pilot local density estimate 
is then used to adapt the amount of smoothing applied to 
each point through ai = ^[fpi^i)/ fp]~^^^ , where fp is the 
geometric mean of the fp{xi). The adaptive kernel estimate 
of the local density fa(x) is then 



27rcr; 



■ exp 



2 



(1) 



The adaptive kernels produce both a higher resolution in the 
clustered regions where it is needed, and increased smooth- 
ing in the low-density regions (see Silverman 1986; Pisani 
1993;1996 for detailed discussion), and the non-parametric 
nature of the method means that the amount of smoothing is 
dependent solely upon the data points rather than any prior 
estimate of what the cluster width should be. The probabil- 
ity density estimate of red galaxies is shown in Fig. |^ with 
the two density peaks, marked A and B, which are assumed 
to mark the two cluster centres. 

Having produced the probability density estimate, each 
galaxy is assigned membership to one of the clusters by fol- 
lowing a path from the original position Xi along the maxi- 
mum gradient of fa{x) until it reaches the local maximum A 
or B. The contribution to the local density from each clus- 
ter, /a (a;), fsix), is then taken to be the sum of the kernels 
for galaxies assigned to that cluster. 

The expected background contamination from red field 
galaxies, fo{x), is unknown, and so we have estimated 
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ID 


Ktot 


A/B 


P{i e 0) 


P(j e A) 


P{i G B) 


1 


15.75 


A 


0.0827 


0.9075 


0098 


2 


16.97 


A 


0.0324 


0.9259 


0.0417 


3 


17.05 


A 


0259 


0.9626 


0.0115 


4 


17.24 


B 


0562 


0.1172 


8266 


5 


17.52 


A 


0.0391 


0.9468 


0.0142 


g 


17.52 


B 


0.0438 


0.0128 


0.9434 


7 


17.60 


B 


0.0319 


0.0159 


0.9522 


8 


18.01 


B 


0306 


0.0164 


0.9529 


9 


18.09 


B 


0.0357 


0.0298 


0.9345 


10 


18.12 


A 


0583 


0.6754 


2663 


11 


18.35 


A 


0.0188 


0.9582 


0.0230 


12 


18.50 


A 


0.6440 


0.2716 


0.0844 


13 


18.59 


A 


0.0200 


0.9595 


0.0205 


14 


18.95 


A 


0.0191 


0.9641 


0.0167 


15 


18.97 


B 


0.0732 


0.2129 


0.7139 


16 


19.09 


B 


0.1860 


0.2110 


0.6030 


17 


19.24 


A 


0.2852 


0.3465 


0.3683 


18 


19.35 


B 


0.2512 


0.1011 


0.6477 



Table 3. Results from the cluster membership analysis of the red 
galaxies. Column A/B indicates whether a galaxy was assigned to 
cluster A or B, and the last three columns show the probabilities 
that a galaxy is isolated or a member of cluster A or B. 

it to be half that of the I — K > 3.75 field galax- 
ies or 0.5 galaxies per square arcminute, as about half 
our I — K > 3.75 galaxies are also classified as red 
(J - A' > 3.75, - / > 2.00). The probability that each 
galaxy is isolated, and hence due to the background com- 
ponent, is just the fractional contribution of the back- 
ground density to the local density P{i G 0) = fo/faixi). 
The probability that galaxy i is a member of cluster A 
is then the fractional contribution of cluster A to the lo- 
cal density, after considering the background contribution 
P(i € A) = {1- P{i e 0))fA{xi)/fa{xi), and for each of the 
red galaxies these probabilities are shown in Table ^. Look- 
ing at this table and Fig. ^ it appears that the two clusters 
are well defined, with 11 out of the brightest 14 galaxies hav- 
ing probabilities greater than 0.90 of having been assigned 
correctly to cluster A or B. Only galaxies #4 and #10 have 
significant doubt over cluster assignment, and #12 appears 
to be isolated, whilst the four faintest red galaxies show no 
clear membership of either cluster, and may only be part of 
the combined structure. 

As a second estimate of the significance of the substruc- 
ture, we examine the null hypothesis that the galaxies are 
all members of a single cluster. One method for examin- 
ing the null hypothesis is the likelihood ratio test statistic 
(LRTS) (Ashman, Bird & Zeff 1994; Kreissler & Beers 1997) 
which evaluates the improvement in fitting the data of a two- 
component model over a single elliptical Gaussian probabil- 
ity density function. The best-fit single elliptical Gaussian 
/(I) (a;) to the data is found using the mean, x, and covari- 
ance matrix. A, of the data, and this unimodal probability 
density function is used to create 10000 bootstrap catalogues 
of 18 galaxies each. The adaptive kernel estimator is applied 
to each bootstrap catalogue in the same way as the origi- 
nal data, and for those catalogues where bimodality is ob- 
served, the cluster membership probabilities P{i e /i), for 
i = 1, N\ fj,= 1,2 are calculated. The probable number of 
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Figure 7. Contour plot of the estimated density distribution of 
optically-red {V - /^2.25, / < 23) galaxies (including all C-M 
relation members) in a 7 X 7arcmin^ (corresponding to 4.06 X 
4.06 Mpc^ at z=1.226) field centred on the quasar. The first 
contour corresponds to a density of 1 galaxy per square arcminute, 
and the separation of successive contours is also 1 galaxy per 
square arcminute. The quasar is marked as previously, and the 
central box corresponds to the field of the K image. 

galaxies in cluster ^ is then given by — y^^., P{i £ /i). 
The best-fit double Gaussians, /(g.,j)(a;), can then be found 
using the means, x^, and covariance matrices, A^, of the 
dataset weighted by P{i € fi) for ^ = 1,2, before normal- 
izing to /(g.jj) (a;)da; = n^. The measure of the fit by the 
g-component model is evaluated through the likelihood value 

jv / a X i-P(»eo) 

ic(5) = n E-As;m)(^0 • (2) 

1=1 \m=1 / 

The evaluation of the improvement in going from a sin- 
gle to a double Gaussian fit is then given by the LRTS 
A = — 21n{I/c(2)/I/c(l)}- The significance of the substruc- 
ture observed is estimated by the probability that a null hy- 
pothesis bootstrap catalogue produces a value of A greater 
than the observed value. For the observed distribution, us- 
ing 10000 bootstrap catalogues, we find that the substruc- 
ture is marginally inconsistent with the null hypothesis at 
the 8.88% significance level. The low level of significance 
is due to the small number of galaxies involved, as many 
of the more significant substructures among the bootstrap 
catalogues were due to just 3 or 4 points that were within a 
few arcsec of one another. 



Clustering across Large-Scales 

As both groups are located near the edges of the K image, 
and so may sufi^er from truncation, the full extent of cluster- 
ing associated with the quasar is estimated by considering 
a 7 X 7 arcmin^ field centred on the quasar from the BTC V 
and I images. Fig. [^indicates that any z>l quiescent galax- 
ies should be amongst the reddest galaxies in V ~ I, and so 
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by selecting those galaxies with V — I > 2.25 the density con- 
trast due to any clustering at z ~ 1.2 should be maximised. 
Note that this selection criterion is slightly redder than the 

V — I > 2.00 selection used previously, in order to reduce the 
contamination of intermediate-redshift quiescent galaxies. In 
the region covered by the K image, only the red galaxies are 
selected, including the 3 with 2.00 < V" - / < 2.25. The esti- 
mated density distribution of V — I > 2.25, 7 < 23 galaxies 
is shown in Fig. |^ and it is clear that the clustering extends 
well beyond the K image. Cluster A now appears to be cen- 
tred on the northern edge of the K image, and cluster B 
appears part of an elongated structure which extends 2-3 
arcmin to the north-east. Two further groups, each of 6-10 
optically-red galaxies, are apparent 3 arcmin south of the 
quasar and 1.5 arcmin apart. The significance of the sub- 
structure formed by clusters A and B appears to be much 
higher now, as the density peaks are moved further apart. 
A re-examination of the null hypothesis, after including the 

V — I > 2.25 galaxies outside the K image, finds the sub- 
structure to be inconsistent with the null hypothesis at the 
1.86% significance level. The increase in significance sug- 
gests that the clusters were truncated by the boundaries of 
the K image, but the increase is also due partly to the larger 
galaxy sample. 



Companion Galaxies to the Quasar Host 

There are two galaxies (labelled CI and C2 in Fig. ^ 
which are only 3 arcsec (12.9ft~^kpc) from the quasar, and 
could well be companion galaxies to the quasar galaxy host. 
Compact companions are found for a significant fraction 
of quasars (e.g. Bahcall et al. 1997) and spectroscopic ob- 
servations confirm that many have stellar populations and 
redshifts within 500 km s~^ of the quasar (e.g. Stockton 
1982; Canalizo & Stockton 1997). It has been suggested (e.g. 
Stockton 1982; Bekki 1999) that these companion objects 
are tidally-stripped cores from galaxies that have recently 
interacted with the quasar host galaxies, and that this in- 
teraction provides an efficient fuelling mechanism for quasar 
activities. Spectroscopic analysis of the companion to quasar 
PG 1700-1-518 (Canalizo & Stockton 1997) finds evidence 
for both a starburst event that occurred roughly 100 Myr 
ago (and so could be coincident with the quasar activation), 
and a relatively old stellar population likely to be from the 
merger progenitor disk. The galaxies CI and C2 are bluer 
than the red sequence galaxies (cf. C3), and so, if associated 
with the quasar, they have undergone a recent episode of 
star-formation, presumably caused by the merger process. 




Separation (arcsec) 

Figure 8. The angular cross-correlation, (9), between red 
{V - I > 2.00,7 - K > 3.75) and blue {V - I < 1.00) galax- 
ies in the field covered by the K image. The errors are assumed 
to be Poissonian with variance Nru- 



To examine the relative spatial distribution of red 
and blue galaxies, the angular cross-correlation function, 
^rb{d) = Nrb{0)/NrR — 1, is determined (see Fig. ^), where 
Nrb{0) is the number of red-blue galaxy pairs with separa- 
tion 9, and NrR{0) is the expected number of pairs with one 
member from the red catalogue and the other from one of 
10000 randomly distributed catalogues. 

At small separations {9 < 20") the red and blue galaxies 
are strongly anti-correlated at the 3o" level, confirming the 
apparent avoidance of the red groups by the blue galaxies in 
Fig. ^. In contrast, for 20" < 9 < 60", there is a correlation 
at the 3(7 level between red and blue galaxies, which is due 
to the 'band' of blue galaxies that bisects the two groups of 
red galaxies. 

The blue galaxies do not contribute to the excess of 
galaxies observed at K < 19, as only 9 of 51 are detected 
in K, the brightest having K = 18.87. These 9 all have the 
red I — K colours indicative of z>l galaxies (e.g. Fig. ^), 
of which 6 help make up the concentration ^30 arcsec from 
the quasar, which suggests that this at least is real. However, 
given that the V — I colour is affected much more by recent 
star-formation than redshift, we cannot say anything about 
the likely redshift of the remaining blue galaxies. 



Blue Galaxies 

To compare this quasar field with the results of HC J we have 
also examined the distribution of blue {V — I < 1) galaxies 
for 7 < 25 (indicated by star symbols in Fig. ^). We find a 
concentration of blue galaxies within 30 arcsec (290/i~^kpc) 
of the quasar, which appears to be extended towards the 
north-east, forming a 'band' that bisects the two groups of 
red galaxies. The band presumably corresponds to a region 
of enhanced star-formation. Few, however, are found near 
the centres of either group of red galaxies. In comparison 
with adjacent fields (over a 7' x 7' region) in the optical 
images we do not find an excess of blue galaxies. 



4 DISCUSSION AND CONCLUSIONS 

We have found a 3.5a excess of 7f < 19 galaxies in the 
2.25 X 2.25 arcmin^ field around the z=1.226 radio-quiet 
quasar 104420.8-1-055739 from the Clowes-Campusano LQG. 
The excess is due entirely to a factor ~11 overdensity of red 
(I — K > 3.75) galaxies, which must have z>0.8 to explain 
their colour. In particular, we find cluster red sequences of 
15-18 galaxies, in the I — K and V — K against K C-M 
diagrams at 7 — Ti" ~ 4.3, V — K 6.9, comparable in both 
colour and magnitudes to red sequences observed for other 
z ~ 1.2 clusters. These red sequences suggest a popula- 
tion of massive ellipticals at the quasar redshift. In the area 
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of the K image we find tentative evidence of substructure 
amongst these galaxies, with two apparent groups 40 arcsec 
to the north and 60 arcsec to the south-east of the quasar. 
An examination of the optical images over a 7 x 7 arcmin^ 
area indicates that this substructure is significant at the 2% 
level, a nd that the clustering extends well beyond the K im- 



ously in a large fraction of cluster galaxies, and they could 
be a major contributor to the Butcher-Oemler effect. Several 
mechanisms have been suggested that could cause firstly the 
triggering and then the termination of a secondary burst of 
star-formation in a galaxy, as a subcluster passes th rough a 



clust er. These include ram pressure from the ICM ( Evrard 



age, forming a large-scale structure 2-3/t ^Mpc across. The 1991), shocks due to collisions between the two ICMs (Eloet- 



overall structure is suggestive of being the early stages of 



formation of a cluster from the progressive coalescence of 
subclusters. 

Only half of the excess oi I — K > 3.75 galaxies is 
accounted for by the red sequence members, and we find 
comparable numbers (15-20) of 'red outlier' galaxies with 
both I - K > 3.75 and V - I < 2.00, which appear to 
fit neither the passively-evolving nor exponentially-decaying 
SFR galaxy models (Fig. Although, all we can say is that 
they are likely to be dusty star-forming galaxies at l<jZ<,2, 
given that such galaxies are found around other z ~ 1.2 
clusters (e.g. TOO; Kajisawa et al. 1999; 2000), and given 
their comparative rarity in field regions, it seems reasonable 
to assume that they are associated with the cluster. This 
would suggest that the Butcher-Oemler effect observed in 
intermediate-redshift clusters increases in strength to higher 
redshif ts, with >50% of likely cluster members exhibiting 



tiger. Burns fc Loken 199(:), and the effect of close galaxy 
encounters ( [Moore'er^l. 1996 ). As the 'band' of blue galax- 



the blu i colours of recent star-formation. 



We also find a concentration of blue {V — I < 1) galax- 
ies within 30 arcsec (130/i~^kpc) of the quasar, with many 
having the red I — K colours of 2^1 galaxies, which suggests 
that the quasar lies in a region of enhanced star-formation, 
in agreement with the results of HCJ. This concentration ap- 
pears extended in such a way as to separate the two groups 
of red galaxies, and it is also notable how the blue galaxies 
appear to avoid the centres of red galaxy clustering. Some 
foreground contamination is likely, but given that such con- 
centrations appear common around LQG quasars at this 
redshift, and the highly significant spatial interrelation be- 
tween the blue and red galaxies, then we can be reasonably 
confident that this concentration is real, and is associated 
with the quasar. However spectroscopic observations of these 
and the other cluster candidates will be required to confirm 
cluster membership, and to provide more quantitative in- 
formation about their star-formation histories, such as their 
approximate ages and current star-formation rates. 

A Possible Cluster Merging Event 

The relative distribution of red and blue galaxies can be 
explained if what we are witnessing is the early-stages of 
merger of the two clusters of red galaxies, which has trig- 



gered 1: oth the band of enhanced star-formation and the 



quasar itself. A compar able distribution has been observed 
for the Coma cluster (Caldwell et al. 199S) with a band 
of post-starburst or 'E-l-A' galaxies located between the 
Coma clust er centre and a se condary X-ray peak. Dynami- 
cal studies ( Burns et al. 195)^ ) indicated that these galaxies 
had passed through the centre of the Coma cluster about 
2 Gyr ago, coincident with the epoch of starbursting pre- 
dicted from the spectra of the post-starburst galaxies. Un- 
usually high blue galaxy fractions have been observed for a 
number of low-reds hift clusters with bimodal X-ray surface 
brightness profiles ( Metevier, Romer fc Ulmer 200C ) , imply- 
ing that cluster mergers can induce starbursts simultane- 



ies is likely to be undergoing or has recently undergone star- 
formation, and because the two groups are relatively close 
together, we suggest that this system is being observed at 
an earlier epoch of the cluster merger process than Caldwell 
et al. (1993), either just before or just after core passage, 
and that the star-formation has been triggered by the in- 
teraction of galaxies with the shock fronts produced by the 
collision of the two ICMs (see Roettiger et al. 1996). 

Cluster merging events are predicted to be relatively 
common at high redshifts in hierarchical clustering 

models (e.g. Press & Schechter 1974; Bahcall, Fan & Cen 
1997; Percival & Miller 1999). Examples of possible merg- 
ing clusters at high redshifts are the CL0023-I-00423 groups 
at z=0.8274 and 0.8452 which, accordi ng to a dynamical 
study, hav e a 20% chance of merging (Lubin, Postman & 



Oke 200C), and the R99 and S97 clusters which are sepa- 



rated by only 2.5/i ^Mpc. 



Comparison with Other Work 

In a study of 7 radio-loud quasars at 1.0 < z < 1.6, Sanchez 
& Gonzalez-Serrano (1999) find excesses of faint [B > 22.5 
and R > 22.0) galaxies on scales of r < 170 arcsec and 
r < 35 arcsec around the quasars, whose numbers, mag- 
nitudes and angular extensions are compatible with being 
clusters of galaxies at the quasar redshifts. In particular 
however, they find that the quasars are in general not lo- 
cated at the peaks of the density distribution, but are some 
40-100 arcsec from them, located on the cluster peripheries. 
This result, and our own, is understandable in the frame- 
work of quasar activity being triggered by the infall of gas 
onto a seed black hole. Firstly, galaxies in the centres of clus- 
ters have previously lost most or all of their gas, by having 
had it stripped off by ram-pressure from the intra-cluster 
medium or by tidal forces from close encounters with other 
galaxies, or by consuming the gas in a starburst during its 
first infall into the cluster. Secondly, the encounter velocities 
for galaxies in the cluster cores greatly exceed the internal 
velocity distribution of the galaxies, making gala xy mergers 
much less effective at triggering nuclear activity (Aarseth & 



Fall 198C) 



In a comparable study of the galaxy environment of the 
radio-loud quasar 1335.8-1-2834 from the z ~ 1.1 Cramp- 
ton et al. (1989) LQG, Tanaka et al. (2000a) obtain re- 
sults which have several similarities to our own. This quasar 
had been part of the HCJ study, and was known to lie in 
a band of blue and emission-line galaxies. Using deep R, I 
and K observations, a number of extremely-red objects with 
the colours of passively-evolving galaxies at the quasar red- 
shift were found, forming a cluster which lies to one side 
of both the quasar and the band of blue and emission-line 
galaxies. They also find a similar population of 'red outliers' 
and estimate the blue galaxy fraction as 60-80%. They also 
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find an indication that this cluster is part of a larger struc- 
ture with groupings of optically-red (i? — / > 1.3) galaxies, 
similar to the clustering near the quasar, found across the 
8 X Sarcmin^ (4.7 x 4.7 h'^ Mpc^ at z=1.086) R,I optical 
images. 

In a wide-field (48 x 9 arcmin'^) optical imaging survey of 
the 1338-1-27 field containing 5 quasars from the Crampton 
et al. (1989) LQG, Tanaka et al. (2000b) detect significant 
clustering of faint red galaxies with 7 > 21 and R — I > 1.2. 
These galaxies are concentrated in 4-5 clusters forming a 
linear structure of extent ~ 10/i~^Mpc that is traced by the 
group of quasars, although only the one radio-loud quasar 
of TOOa appears to be directly associated with any of the 
rich clusters. The immediate environments of the other four 
LQG quasars, all radio-quiet, appear relatively poor in terms 
of red galaxies, although three have excesses of blue and 
emission-line galaxies (HCJ95) indicating that they are lo- 
cated within regions of enhanced star-formation. This is the 
clearest evidence yet that LQGs trace large-scale structure, 
even if the majority of the member quasars are only directly 
associated with regions of enhanced star-formation, rather 
than rich clusters. 



Interpretation - Mechanisms for Quasar Formation 

As several of the HCJ quasars were found in regions of en- 
hanced i Btar formation, and ac both tho quaear of TOOa and 



this pajkpT- ar-P InfatpH in 'hanrlg' nf pnhanfprl gtar-frirmatinn 

in between, or on the peripheries of, clusters, we propose a 
causal link between the quasar and star-formation whereby 
both are triggered by the same mechanism; the interaction 
between the galaxy and the intra-cluster medium. If a galaxy 
can be disrupted sufficiently by its passage through the ICM 
to cause it to undergo starbursting, then if it also contains 
a supermassive black hole, enough gas may be channelled 
onto the nucleus to trigger a phase of quasar activity. Such 
a mechanism explains the finding of quasars in regions of 
enhanced star-formation more naturally than the galaxy 
merger model, as it allows many galaxies to be affected si- 
multaneously, although there is good evidence that a large 
fraction of quasars have been triggered by galaxy mergers. If 
there is a connection between the quasar and star-formation 
activation, then it is likely that these quasars are found pref- 
erentially in clusters with high blue galaxy fractions, and 
that these are not representative of « ~ 1.2 clusters as a 
whole. 

This work and previous studies show that searching for 
sources with the / — A" > 3.75 colours characteristic of qui- 
escent galaxies at z^l is an efficient means of locating 
clusters. By adding a second optical band it is then possible 
to obtain qualitative information on the star-formation his- 
tory of these galaxies. Using this information along with the 
relative spatial distribution of quiescent and star-forming 
galaxies, a more complex picture arises in which the evolu- 
tion of galaxies, quasars and clusters are all interrelated. 
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